Volume 78, Number 17 



PHYSICAL REVIEW LETTERS 



28 April 1997 



Anomalous Size Dependence of Relaxational Processes 
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We consider relaxation processes that exhibit a stretched exponential behavior. We find that in those 
systems, where the relaxation arises from two competing exponential processes, the size of the system 
may play a dominant role. Above a crossover time f x that depends logarithmically on the size of the 
system, the relaxation changes from a stretched exponential to a simple exponential decay, where the 
decay rate also depends logarithmically on the size of the system. This result is relevant to large-scale 
Monte Carlo simulations and should be amenable to experimental verification in low-dimensional and 
mesoscopic systems. [S003 1-9007(97)02972-4] 



PACS numbers: 64.60.-i, 64.70.Pf 

In recent years it has become clear that many relaxa- 
tional processes in macroscopic systems can be character- 
ized by a relaxation function Q(t) that exhibits a stretched 
exponential behavior 



e(0~G(0)exp[-(f/ry], 



(1) 



where < (3 < 1 . Examples include viscoelastic relaxa- 
tion [1], dielectric relaxation [2], glassy relaxations [3- 
5], relaxation in polymers [6,7], and long-time decay in 
trapping processes [8]. Many more examples [9-13] 
suggest that Eq. (1) is common to a very wide range of 
phenomena and macroscopic materials. 

The origin of the stretched exponential is not always 
clear. In many systems it is assumed to be the result 
of a competition between two exponential processes. In 
some cases, e.g., trapping processes at long times, this 
assumption is well established, while in others, such 
as relaxation in glassy materials, this assumption has 
been controversially discussed [14,15] and alternative 
models have been also suggested [10,16-18]. Less is 
known, both experimentally and theoretically, on the 
corresponding behavior in mesoscopic systems where we 
expect the relaxation to depend on the system size. 

In this Letter we argue that if the stretched exponen- 
tial is due to two competing exponential processes, there 
exists a characteristic time r x , which depends logarithmi- 
cally on the size of the system, above which there is a 
crossover to an exponential decay. Thus, by varying the 
size of the system this crossover time changes. This can 
serve as an experimental test for identifying the origin of 
the mechanism leading to stretched exponential decay. 

We assume that the relaxation function of the whole 
system can be represented by an integration over all 
possible states n, namely, 



Q(t) 



Jo 



<$>{n)Q{n, t) dn . 



(2) 



Here, O(n) is the probability that state n is occupied and 
Q(n,t) is the dynamic relaxation of the nth state. 



Usually, in the case of a stretched exponential behavior, 
<J>(n) is assumed to behave as O(n) ~ exp( — an a ), while 
Q(n,t) decays exponentially with time as Q(n,t) ~ 
exTp( — bt/n y ). A number of dynamical models that yield 
a stretched exponential decay can be formulated in terms 
of Eq. (2). These include the long-time behavior in the 
trapping problem [8], the target problem [18], direct 
energy transfer [18], hierarchically constrained dynamics 
[14] and others. We now concentrate on two examples: 
The first example is a particle diffusing in a J-dimensional 
system with randomly distributed traps, where we are 
interested in the survival probability Q(t) of a particle. 
Here the state n represents a particle in a trap-free region 
of linear size n; <&(«) is the probability for the occurrence 
of a size n trap-free region, and Q(n,t) is the survival 
probability of the particle in this region [8]. The exponent 
a is the dimension d of the system, and y = 2 due to the 
diffusional motion. The second example is hierarchically 
constrained dynamics, a model that has been proposed to 
account for glassy relaxation [14]. This model assumes 
that the relaxation of level n populated by spins occurs 
in stages, and the constraint imposed by a faster degree of 
freedom must relax before a slower degree of freedom can 
relax. This implies that the time scale of relaxation in one 
level is subordinated to the relaxation below. A possible 
realization considered in [14] and here is a system with a 
discrete series of levels where the relaxation time of level 
n is t„ ~ n y [corresponding to the exponential form of 
Q(n,t) in Eq. (2)], and the weight factor of level n is 
O(n) ~ e~ an [12], corresponding to a = 1. The first 
exponential in Eq. (2) is, accordingly, the probability to 
occupy level n and the second exponential represents the 
decay of that level. 

We can evaluate the long-time behavior of the integral 
in Eq. (2) using the method of steepest descent. The 
main contribution to the integral arises from the maximum 
of the integrand in (2), which is obtained from the 
minimum of the function — an" — bt/n y appearing in 
the exponent. This yields that the main contribution to 
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(2) comes from 



n* = {ybt/aa) l/(a + y \ 



(3) 



leading to Eq. (1) with f3 = a/(a + y) < 1, and r = 
(a/by)a-y/ a [y/(y + a)] 1+y / a . 

However, as we show below, these arguments are valid 
only in the thermodynamic limit where the system size 
is infinite. For a finite system with a finite number N 
of spins (in the hierarchical constraint system) or a finite 
number N of traps (in the trapping system), the relaxation 
function depends explicitly on N. Since our discussion 
is quite general for systems described by Eq. (2), in what 
follows we refer to spins and traps in the above examples 
as elements. 

For a single finite system consisting of N elements, the 
relaxation function Q{t) represents an average quantity 
over the N elements 



Q{t) = - J«(»)fi(B,/), 

W 



(4) 



where the sum is over all possible states n and m(n) is 
the number of elements at state n, with Y.{ n } m ( n ) = N. 
Since the sum in (4) is over exponential functions, the 
value of Q(t) will fluctuate for different sets of N. There 
will be a distribution of Q(t), and we are interested in the 
typical Q(t), which is around the peak of this distribution. 

In the thermodynamic limit N — ► °°, all states n are 
occupied, m(n)/N can be identified with and 
Eq. (2) follows. For N finite, in contrast, there exists 
a characteristic "maximum" state n = n max (N), and this 
"max should replace the upper limit (°°) in Eq. (2), 



Q(t) 
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<£>(n)Q(n,t) dn . 



(5) 



To estimate how n max depends onAf, we note that the 
typical number of states n in a sample of N elements is 
Z(n) = N<Z>(n) = Nexp(-an a ). States with Z(n) « 1 
will not occur in a typical system of N elements, and this 
yields 



lnN\ l/a 



(?) 



(6) 



If <SC re max , the upper limit in (2) can be approxi- 
mated by infinity and thus leads to Eq. (1). However, if 
ri* » ra max the main contribution to Eq. (5) will not be 
from the maximum of the integrand, which is outside the 
range of integration, but from n max . Thus, for n* » w max 
we expect 



Q(t) = Q(0)e- ht ^\ 



(7) 



where the time constant of the relaxation, scales as 

(lnN) y / a . The crossover time from a stretched exponen- 



tial [Eq. (1)] to an exponential [Eq. (7)] can be estimated 
from the condition n* = « max , from which follows 



aa /InN 
yb V a 



\ + y/a 



(8) 



The striking point in Eq. (8) is the logarithmic depen- 
dence on N, which puts t x in the range of observable 
time scales measurable in mesoscopic systems. Indeed, 
the corresponding relaxation value Q (tx) scales as 



Q(t x ) ~ N~ a /y, 



(9) 



independent of the microscopic parameters a and b. In 
the case of the trapping relaxation mechanism where 
a = d and y = 2 we obtain 



Q(tx)/Q(0) ~N- d '\ 



(10) 



Q(tx)/Q(0) ~ N d ^ 2 , while in the hierarchical constraint 
dynamics 



G(fx)/G(0) ~ AT 



-1/7 



(ID 



It is known [8(e)] that in both examples, for an infinite 
system, the stretched exponential behavior of Eq. (1) sets 
in only at very long times. Thus we expect that in 
the finite system, the crossover will mask the stretched- 
exponential pattern. 

To test our analytical approach, we performed Monte 
Carlo simulations on both the trapping model and the 
hierarchical constraint model. In the trapping model, 
we consider one- and two-dimensional systems with a 
fixed concentration c = 0.5 of randomly distributed traps 
and vary the size N/c of the system. We calculated 
numerically the survival probability Q(t) of a particle 
as a function of t and N. In the hierarchical model 
we have chosen t„ ~ n, i.e., y = 1. We calculated 
the relaxation function for system sizes varying from 
N = 10 2 to Af = 10 5 . 

As mentioned earlier, the relaxation function fluctuates 
for different sets of N. For obtaining the typical 
behavior of Q{t), we have considered therefore the 
"typical" average <2(?) ty p = exp((ln Q(t))), where the 
brackets denote an average over many sets of N elements 
[19]. For simplicity, we shall drop the index "typ" in the 
following. 

Figure 1 shows — ln[<2(?)/<2(0)] as a function of t in a 
double logarithmic plot for (a) the trapping model in d = 
1 and d = 2, and (b) the hierarchical constraint model, 
both for several system sizes. In all cases, a crossover 
from an exponent f3 < 1 (at small t) towards /3 = 1 (at 
large t) can be easily recognized. The crossover time tx 
shifts towards larger values when N increases. 

To study the crossover behavior in a more quantitative 
manner, we have plotted in Fig. 2 the local exponents (3 
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FIG. 1. Plot of -ln[62(0/G(0)] as a function of t in a 
double logarithmic presentation for (a) the trapping model in 
d = 1 and d = 2, and (b) the hierarchical constraint model, 
for several system sizes. For the trapping model, the system 
sizes are N = 2 X 10 3 (open square), 2 X 10 5 (open circle), 
2 X 10 7 (open up triangle), 2 X 10' (open down triangle) in 
d = 1, and N = 9 X 10 2 (full square), 9 X 10 4 (full circle), 
9 X 10 6 (full up triangle) in d = 2. For the hierarchical model, 
the system sizes axe N = 10 2 (full square), 10 3 (full circle), 10 4 
(full up triangle), 10 5 (full down triangle). 



obtained from the local slopes of Fig. 1, as a function of 
t. In both systems, for a fixed system size N, (3 first 
decreases with t, reaches a minimum value at a certain 
time that can be identified with f x , and then increases 
monotonically with time towards (3 = 1. The figure 
shows that the minimum value of (3 has not yet reached 
its asymptotic value predicted for infinite systems, i.e., 

= 1/3 (d = 1) and /3 = 1/2 (d = 2) for the trapping 
system and (3 = 1/2 for the hierarchical system. 

To show the dependence of the crossover time t x on 
the system size A' we have plotted, in Fig. 3, the values of 

t x as a function of In A'. The crossover time was 

obtained numerically from the position of the minima of 
the curves in Fig. 2. The resulting straight lines are in full 
agreement with the prediction of Eq. (8), supporting our 
analytical approach. 
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FIG. 2. Plot of the local exponents (3 calculated from the 
successive slopes of the corresponding curves in Figs. 1(a), for 
the trapping model and 1(b) for the hierarchical model. The 
horizontal dashed lines represent the corresponding asymtotic 
(jV — » oo, t — ► oo) values of f3. 



In the following we discuss the relevance of our results 
to Monte Carlo simulations and experiments. There exists 
a long standing puzzle in Monte Carlo simulations of 
the trapping problem in d = 2 and 3, that the predicted 
stretched exponential could not be observed [8], even for 
survival probabilities Q(t)/Q(0) down to 10" 21 in d = 2 
[8(b)] and 10" 67 in d = 3 [8(g)]. 

Our finding of the logarithmic dependence of Q{t) on 
the system size N explains this puzzle. The Monte Carlo 
simulations in d = 2 and 3 were typically performed 
on 10 3 configurations with about 10 4 traps, which is 
equivalent to having a single system with A' ~ 10 7 
traps. Using Eq. (10), we expect for N = 10 7 traps 
Q(tx)/Q(0) = 10" 7 in d = 2. Indeed, for times above 
tx the exponent (3 approaches unity as predicted by 
our theory and as seen clearly in Fig. 2(a). Moreover, 
for this system size f3 never reaches the predicted 
thermodynamic value (3 = 0.5, the minimum value of 
13 is about 0.65. For d = 3, Q(t x )/Q(0) = 10" 11 thus 
for smaller survival values (t > t x ) one again expects 
increasing values of (3 approaching unity. This explains 
the exponential decay found in the early Monte Carlo 
simulations. Our results show that this is not an artifact 
but due to the finite size of the system. Moreover, they 
clearly indicate that the thermodynamic limit cannot even 
be reached in one-dimensional macroscopic systems. 

It would be of interest to test the above prediction experi- 
mentally by preparing experimental realizations where size 
effects can be controlled. Equations (8) and (10) suggest 
that the behavior around the crossover can be measured ex- 
perimentally. For the trapping problem in linear systems, 
which has been studied experimentally [20,21], we expect 
for 10 8 sites and concentrations of traps c between 10~ 4 
and 10" 2 , that g(f x )/G(0) ~ 10" 2 -10" 3 , which is a sur- 
vival range that can be detected experimentally. The same 
arguments are valid for the target problem and therefore a 
similar crossover from stretched exponential to exponen- 
tial decay is expected in relaxation experiments in low- 
dimensional geometries [22]. Mesoscopic systems such 
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FIG. 3. Plot of tx (a+y) as a function of lnN, for (a) the 
trapping model and (b) the hierarchical model. The straight 
line supports Eq. (8). The crossover times t x were obtained 
from the positions of the minima of Fig. 2. 
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as quantum dots, are also promising candidates for ex- 
periments where the crossover can be relevant. Identify- 
ing the logarithmic size dependence in experiments may 
provide support to the theories claiming that the observed 
stretched exponential is due to competing exponential pro- 
cesses, represented by Eq. (2). 

We would like to thank Julia Drager for valuable 
discussions. This work was supported by the German 
Israeli Foundation (GIF). 
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